Comment on "Canonical and Mircocanonical Cal- 
culations for Fermi Systems" 

In the context of nuclear physics Pratt jij recently in- 
vestigated nonintcracting Fermi systems described by the 
microcanonical and canonical ensemble. As will be shown 
his discussion of the model of equally spaced levels con- 
tains a flaw and a statement which is at least confusing. 

In his Letter Pratt gives an alternative derivation of 
an expression (Eq. (3) of Ref. M) in which the canonical 
partition function Za{P) of a system with A fermions at 
inverse temperature (3 is expressed in terms of the parti- 
tion function with fewer particles Z J 4_„(/3) Starting 
with Zq((3) = 1, Za(/3) can thus be calculated recur- 
sively. This result also leads to an expression in which the 
canonical occupation probability (rik} c for a one-particle 
state with energy is expressed as a sum over terms 
which contain Za-u(P) and exp(— n[3ek) Za(0) 
can be expressed as the (discrete) Laplace transform of 
the number of states NA,E r of a system with A particles 
and fixed energy E r . Here the index r indicates, that for 
a finite system the energy is discrete and given as the 
sum over A one-particle energies Furthermore the 
product Za(P) (nk) c is given as the Laplace transform 
of Na.e t i n k) mc where (nk) mc is the microcanonical dis- 
tribution function. Starting with the expressions for the 
canonical ensemble and using inverse Laplace transfor- 
mation it is thus possible to obtain recursive relations 
for Na,e t Band ( n k) mc - They are given in Eqs. (7) and 
(8) of Ref. II . 

As an example Pratt investigated a model with equally 
spaced (spacing g) nondegenerate fermionic levels. The 
A-particle groundstate energy is given by Eq(A) — 
gA(A + l)/2 and the total energy E r (A) can be writ- 
ten as E r (A) = E Q (A) + gr, with r = 0, 1,2, . . .. His 
discussion suggests that there is a direct connection be- 
tween the number of states with an excitation energy gr 
and the number of partitions of the integer r into dis- 
tinct parts. This is not the case. Instead the model can 
be related to the unrestricted partition problem [Qg, i.e. 
the problem of finding the number of ways the integer 
r can be written as a sum of smaller integers, where it 
is allowed to repeatedly use the same integer. The very 
fact that the same integer can be used more than once is 
essential in connection with the mapping of the fermionic 
system on a bosonic one, a procedure which is known as 
bosonization B. For r < A, NA,E r is given by the num- 
ber of partitions p(r) of r which is independent of A. For 
r > A this mapping breaks down and Na,e t < p{ r ) 
In Fig. 3 of his Letter Pratt presents data for (n^) and 
r = 20 calculated using his Eqs. (7) and (8). He does not 
give the number of particles for which the data have been 
obtained. We have been informed by Pratt that r < A, 
which implies that the data have to be independent of A. 
For r < A the microcanonical occupation probability ful- 
fills the symmetry relation (nk)„ lc + (riA-k+i) mc = 1 and 
is a monotonic function of k Mja . These observations are 


essential in connection with the question how the grand 
canonical distribution function (the Fermi function) is 
appoached in the thermodynamic limit. Both require- 
ments are apparently violated by the data presented by 
Pratt. We have been informed by the author that in 
Fig. 3 he accidentally used the wrong set of data. The 
data shown have been obtained for a model of randomly 
placed levels. Fig. [j] shows the correct microcanonical oc- 
cupation probability for the equally spaced level model, 
r = 20, and different system sizes. To obtain the data for 
r < A we have used the three different algorithms pre- 
sented in Refs. @], ||, and Q which all lead to the same 
curve. For r > A the formulas presented in Refs. Q and 
H cannot be used and only the new result of Pratt which 
does not make direct use of the mapping to the partition 
problem is applicable. {nk) mc displays monotonic behav- 
ior also for r > A. As a test of the validity of our data we 
have always checked the sum rule (nk) mc = A which 
is fulfilled up to very high precision. 
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FIG. 1. The microcanonical occupation probability (Hfc) mc 
for excitation energy 20g and different system sizes. For com- 
parison the Fermi function for temperature r = ksT/g = 3.49 
is shown as the solid line (see Ref. [1]). 
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